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Abstract 

From the expression for the electromagnetic field in the neighborhood of a point charge we 
determine the rate of electromagnetic momentum flow, calculated using the Maxwell stress tensor, 
across a surface surrounding the charge. From that we derive for a ‘point charge’ the radiation 
reaction formula, which turns out to be proportional to the first time-derivative of the acceleration 
of the charge, identical to the expression for the self-force, hitherto obtained in the literature 
from the detailed mutual interaction between constituents of a small charged sphere. We then 
use relativistic transformations to arrive at a generalized formula for radiation reaction for a point 
charge undergoing a relativistic motion. 
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I. INTRODUCTION 


The radiation reaction, first proposed by Lorentz,—— was later derived in detail by 
Abraham^ and Lorentz^ from the self-force of an arbitrarily moving small charged sphere. 
It is in fact remarkable that Abraham’s formula,- derived before the advent of the special 
theory of relativity, includes all the relativistic terms correctly. Von Laue derived the same 
from a relativistic generalization of Lorentz’s non-rclativistic formula for radiation reaction 
by expressing it in a 4-vector form.- An exhaustive treatment of these formulations was 
later provided by Schott^ and Page and Adams.- The computations of the self-force of a 
charged sphere have a long history.-^ These formulas are now available in various forms in 
modern text-books.- 1 ^— In these derivations one usually considers for the charge particle, 
a classical uniform sphere model (a solid ball or a spherical shell) of radius r D , which may 
be made vanishingly small in the limit. One calculates the force on an infinitesimal element 
of the charged sphere due to the fields from all its remainder parts, with the positions, 
velocities and accelerations of the latter calculated at the retarded times. Then the net 
force on the charge is calculated by integration over the whole sphere. To make it simple, 
the force calculations are usually done in the instantaneous rest-frame of the charge where 
it is generally assumed that the motion of the charged varies slowly so that during the light- 
travel time across the particle, any changes in its velocity, acceleration and other higher 
time derivatives are relatively small. Then one can keep only linear terms v/c, vr Q /c 2 , 
vr 2 /c 3 , etc., in the formulation of self-force, which turns out to be proportional to the rate 
of change of acceleration. Further the self-force is found to be independent of the radius of 
the small sphere. An alternative derivation of the radiation reaction equation was given by 
Dirac,— that involved both retarded and advanced solutions, which is sometimes open to 
criticism based on the causality arguments. The problem of radiative losses and its relation 
with radiation reaction on the charge has been extensively discussed in the literature.—"^ 

The question that we examine here is: could the radiation reaction expression, which is 
independent of the radius of the small sphere, be actually derived in the case of a point 
charge? Of course one cannot meaningfully talk of a self-force on a point charge due to 
electromagnetic interaction between its constituents. But there is a possibility to exploit 
the law of conservation of mechanical and electromagnetic momentum of a system comprising 
a charge and its electromagnetic fields. We shall endeavor to do so here. We accordingly 
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show that the expression for radiation reaction derived this way is exactly the same as has 
been derived for a sphere of vanishingly small radius. 

We shall first examine electromagnetic field in the neighborhood of an arbitrarily moving 
point charge. From that we shall calculate the radiation reaction, for a non-stationary charge 
moving with a non-relativistic velocity, from the electromechanical momentum conservation, 
using Maxwell stress tensor. Later we will transform the radiation reaction formula to a 
relativistic case. 


II. ELECTROMAGNETIC FIELD IN THE NEIGHBORHOOD OF A POINT 
CHARGE 


The electromagnetic field (E. B) of an arbitrarily moving point charge e at any time t Q 
is given by—^ 


E 


e(n — v/c) 
7 2 r 2 (l — n • v/c) 3 

B = 


e n x {(n — v/c) x v} 
rc 2 (1 — n • v/c) 3 

n x E, 


( 1 ) 

( 2 ) 


where the square bracket indicates that quantities like velocity (v), acceleration (v), and 
the Lorentz factor (7 = 1 / \/(l — (v/c) 2 )) of the charge are all in terms of their values at the 
retarded time t Q — rjc. The radial vector r = rn is from the retarded position of the charge 
to the field point where electromagnetic fields are being evaluated. 

Using the vector identity v = n(v.n) — n x {n x v}, we can decompose the electric field 
in Eq. (JTj), in terms of the radial (along n) and transverse components as^S^ 


E = 


e n 


7 2 r 2 (l — n • v/c) 2 


enx {(n — v/c) x (yv + y 3 vr/c)} 
7 3 r 2 c (1 - n • v/c) 3 


( 3 ) 


- ret 


the quantities on the right hand side evaluated at the retarded time t Q — r/c. This is a 
general expression for the electric field of a charge, with the radial (along n) and transverse 
(perpendicular to n) terms fully separated. 

Initially we assume that motion of the charge is non-relativistic and that it varies only 
slowly so that any change in the velocity, acceleration etc. during the light-travel time across 
the region of interest are relatively small. This means that |v| -C c, |v|r/c c, |v|r/c 
|v|, etc. Therefore we drop non-linear terms in v, v etc., in our formulation. 
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Thus keeping only linear terms in v and v, from Eq. (j3J) we can write 



[1 2n • v] 


e n 


+ e n x 



ret 


n x (v + vr /c) 


r 2 c 


ret 


( 4 ) 


Now, under the above conditions, one can make a Taylor series expansion of the retarded 
velocity and acceleration vectors around the present time t Q 

„2 


v Q r v Q r 

v = v “- — + i? + 


v Q r 

v = v G -1-, 

c 


(6) 

( 6 ) 


all vector quantities on the right hand side evaluated at the present time t Q . 

Substituting for v and v from Eqs. ([5]) and © into Eq. ((4j) , and dropping terms of order 
r or its higher powers, that will become zero as r —> 0, we get for the electric field in the 
neighborhood of an accelerated point charge as 

‘n X v r , n x vT 


E = e n 


1 2n ■ v n 


r 2 c 


2n • v Q n • v Q 

o ' o 


rc* 


and the magnetic field as 


B = e 


n x v 0 

r 2 c 


+ 


+ e n x 


n x v n 


r 2 c 


2c 3 


2c 3 


( 7 ) 


( 8 ) 


It should be noted that in Eqs. (J7D and ([HD, while the radial distance r and the unit vector n 
are specified with respect to the retarded position of the charge, the velocity and acceleration 
at the retarded time have been expressed in terms of the charge motion at the present time 
t D through Eqs. © and (j6j) . 


III. THE COMPUTATION OF RADIATION REACTION FROM MAXWELL 
STRESS TENSOR 


Law of conservation of momentum in an electromagnetic system, containing both charges 
and fields, states that at any given instant, say t 0 , the rate of change of the combined 
momentum in a volume V, comprising the mechanical momentum of the charges (P mec h) 
and the electromagnetic momentum of the fields (Pfieid), is equal to the rate of momentum 
flow into the volume through a surface £ surrounding that volume.— 1 ^ 


d(Pmech T Pfieid) 
(It 


d£ n • T • 


( 9 ) 
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where all quantities are calculated for the same instant of time t a . The vector component 

■O- -f->- 

n • T of the Maxwell stress tensor T, represents the radial momentum flow into the volume 
per unit area of the surface, and can be calculated fro m 1 °' 13 


n • T 


1 

47T 


(n • E)E + (n • B)B 


n(E 2 + B 2 )l 
2 . 


( 10 ) 


The volume integral of the electromagnetic field momentum enclosed within the spherical 
surface £ is given by 

Pfieid = ^/dl/(ExB). (11) 


Let us consider a charge moving with an arbitrary but non-relativistic speed (v -C c). 
We assume that the moving charge was at time t = 0 at a point O, and now at time t Q 
the charge is at its ‘present’ position P. Let £ be a spherical surface of radius r = ct Q , 
centered on O and let V be the volume enclosed within (Fig. 1). The electromagnetic 
field at any point on surface £ is determined from the position and motion of charge at O, 
which is the corresponding time-retarded position of the charge. We want to use the law 
of momentum conservation to compute the mechanical force on the charge at P. It should 
be point out here that although the present position of the charge is not at the center of 
the sphere £, yet the charge and its neighborhood is still enclosed well within the surface 
£ (as v <C c) and the electromagnetic momentum conservation is as much applicable to the 
electromagnetic momentum flux through this surface as to a surface centered on the charge. 
As the charge moves, different spherical surfaces within the volume are causally related to 
different retarded positions of the charge, and the corresponding radial vector r and the unit 
vector n are constantly changing. In Fig. 1, we have shown another spherical surface £i of 
radius r —dr = c(f 0 —dt), centered at the corresponding time-retarded position of the charge 
at a distance vdf from O. In order to calculate the volume integral of the field momentum, we 
first evaluate the volume element between the two spherical surfaces £ and £i, centered on 
two different retarded positions of the charge, separated by vdt (Fig. 1). The radial distance 
between the two spherical surfaces is not a constant dr but is instead dr(l — n- v/c), with the 
corresponding volume element given as dV = 27rr 2 (l — (v/c) cos 6) sin 6 dr d6. Using Eqs. (171) 
and ()8j) , and dropping all terms which are non-linear in v Q or its derivatives, or those which 
upon integration over the solid angle of the whole sphere yield a zero net value, or those 
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n-T 



FIG. 1. Geometry for calculating the volume element in the case of a moving charge. The spacing 
between the two spheres in direction n is dr{l — n • v/c). The vector component of the Maxwell 
stress tensor n • T represents the radial electromagnetic momentum flow into the volume through 
a unit area of the surface X. 


proportional to r or its higher powers which would vanish when r —>• 0, we get 
— n2 '■ nx (nx v n ) 2e 2 v n r r dr 2e 2 v n fl T 


P fiplrl 


47TC 


dP 


r 4 c 


3 c 2 




3 c 2 


r 


( 12 ) 


This is the electromagnetic held momentum in the volume V, enclosed within the spherical 
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surface £, at time t Q . The lower limit of the integral is restricted to e, as done even in static 
Coulomb fields to avoid divergence of the self-field energy at low r. Then we get 


dPfield _ 2e 2 Vp ri _ r 
d t 3 c 2 _e r. 


(13) 


The first term on the right hand side in Eq. (fT3l) represents the change in momentum of the 
self-fields of the charge as its velocity changes. Its divergence, as e —» 0, is unavoidable due 
to the ‘point’ nature of the considered charge distribution where the Coulomb field energy 
of even a stationary charge diverges. The formulation for a point charge remains arbitrary 
to that extent.— 

We substitute for electric and magnetic fields from Eqs. (J7|) and (JHJ) into Eq. ffTUl) . and 
integrate over the surface £ to obtain the electromagnetic momentum flow into enclosed 
volume V at t Q as 


d£ n • T 


— <f r 2 dft 
4(/r Jn 


n 


2n • v n 2n • v n 


■r 2 c 


rc 


+ 


n • v n 


n 

H rr x 


f n x v 0 

n x v Q V 

l r 2 c 

2c 3 jJ 


(14) 


where is the solid angle. Upon integration, terms containing v Q get cancelled, leaving us 
with 


d£ n • T = — <f dfi 
s 47 t Jn 


2n(n • v 0 ) n(n ■ v Q ) + W 

o ' 


rc* 2c 3 


= — / dO sin 9 

2 Jo 
2e 2 . 2e 2 .. 

3rc 2 ° 3c 3 ° 


2v 0 cos 2 6 v 0 cos 2 6 + v Q 

o + 


rc^ 


2c 3 


(15) 


Now from Eq. (JUJ), the mechanical force F = dP mec h/df, on the charge enclosed within V at 
t 0 is 

F = / s dEn ' (16) 
Substituting Eqs. ()T3j) and f)T5|l into Eq. (ITS]) , we get the electromagnetic force on the charge 
at time t a as 


F = 


2e 2 


2e 2 .. 


(17) 


3 e c 2 ° 3c 3 ° ’ 

where v 0 and v Q respectively are the acceleration and the rate of change of acceleration of 
the charge at t Q . An interesting thing to note here is that one could have arrived at Eq. 
without evaluating the field momentum (Eqs. (fl2|) and ()T3j) ). by using e as the radius of the 
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spherical surface £ directly in Eq. (fT5]l . However, from Eq. (IT2]) we learn that 2e 2 v Q /(3 c 2 e) 
is the held momentum in the volume exterior to a radial distance e from the point charge 
(we let r —* oo in the last term in Eq. (1T2|) ). due to its present velocity v Q . This is true 
for any value of e. This of course is a well-known result for a uniformly moving charge,— >22 
but now we know that this formula is applicable even if the charge were moving arbitrarily 
but with a non-relativistic motion. As the charge velocity changes due to its acceleration, 
the held momentum, oc v G , also must accordingly change. Therefore, from the functional 
form of the hrst term on the right hand side of Eq. ffTTjl . it is evident that this term is due 
to nothing else but the changing electromagnetic momentum in the helds because of the 
changing velocity of the charge and has nothing to do directly with the radiative losses by 
the charge. It is only the second term in the force formula Eq. f|17p . independent of r, which 
represents radiation reaction and is exactly the same as obtained from a derivation of the 
self-force of an accelerated charged sphere of small dimensions.-^-^— 

Thus the expression for radiation reaction 


F = 


2e 2 .. 


(18) 


being proportional to rate of change of acceleration, has been derived here from the Maxwell 
stress tensor of a point charge. 

From the Poynting flux in the neighborhood of a point charge moving with a non- 
relativistic but otherwise arbitrary velocity v Q , it has recently been shown^i that the ra¬ 
diative power loss is given by 



(19) 


which, using Eq. (jTSlh can be written as 


V = -F • Vo . 


( 20 ) 


This is consistent with F being a radiative drag force on the moving charge. 


IV. RADIATION REACTION FOR A RELATIVISTIC MOTION OF THE CHARGE 

We calculated the radiation reaction (Eq. fjl8|l ) for a point charge having a non-relativistic 
but otherwise arbitrary motion. We shall now use relativistic transformations to obtain 
expression for radiation reaction for a charge moving with relativistic speeds. 




In an inertial frame 1C (lab-frame!) let the charge at some instant t be moving with 
velocity v and acceleration v and let at that event JC be the instantaneously rest frame 
of the charge, thus moving with a velocity vy = v with respect to 1C. The proper time 
interval d tf in rest frame 1C' is related to the corresponding time interval in the lab-frame 1C 
as d t = yd tf. 

In the instantaneous rest frame 1C' we rewrite radiation reaction (Eq. ([T8]) ) as 


F' = 



( 21 ) 


where we have dropped the subscript from v D and its temporal derivatives. We now want 
to transform this relation to the frame 1C. 

The force components transform from the rest frame 1C' to the lab-frame 1C as^2 


Fj, = F|| , (22) 

K = yFx . (23) 

The parallel component of acceleration transforms as 

vj, = 7 3 V|| , (24) 

from which we get 

v|| = 7(3 7 5 (v ■ v/c 2 )v|| + 7 3 V||) , (25) 

On the other hand the perpendicular component of acceleration transforms as22 

v(l = 7 2 v ± + y 4 (v ■ v/c 2 )vi . (26) 


From that we get 


Vj_ = 27 5 (v • v/c 2 )vi + 7 3 vx + 7 5 (v • v/c 2 )v ± 


(27) 


It should be noted that though the instantaneous value vj_ = 0 (by the definitions of frames 
1C' and 1C), we retained it in Eq. (126|) as it is needed in the derivation of Eqs. ([27]) . where 
vi need not be zero. 

Thus for the parallel component of force from Eqs. (1TT|) . (1771) and (175]) we get 


Fii = 


2c 2 

3c 3 


37 6 (v-v)v|| 4 „ 

-7S-+ 7 V H 


(28) 
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while for the perpendicular component of force Eqs. (ED, El and (1271) lead us to 


f 

X 3c 3 


37 4 (v-v)v ± 2 .. 

-5-+ 7 v ± 


( 29 ) 


Combining the two force components and after a rearrangement of various terms we get the 
required relativistic expression for the radiation reaction as 


2_,2 


F = 


2 e 7 
3c 3 


7 2 (v.v)v , 3 7 2 (v'v)v . 3 T 4 (v-v) 2 v 


t- <7 c- (30) 

which is the relativistic form of the radiation reaction formula, earlier derived in literature 
by various different methods.— - — 

Frorn the relativistic expression (Eq. (130)1 ) for the radiative drag force, the radiative power 
loss in a relativistic case is 


2-,4 


V = —F • v = 


2 e 7 ' 
3c 3 


2 (v-v) 2 

V ■ V + 37 2 v ’ 


Cf 


(31) 


V. CONCLUSIONS 

From the electromagnetic held in the neighborhood of a ‘point charge’ in arbitrary motion, 
we derived here the radiation reaction from the rate of electromagnetic momentum how, 
calculated using the Maxwell stress tensor, across a surface surrounding the point charge. 
Our result that the radiation reaction is proportional to the hrst time-derivative of the 
acceleration of the charge, is in complete agreement with the self-force of a small charged 
sphere, determined in literature from a detailed mutual interaction between its constituents. 
We also derived a generalized formula for radiation reaction for a relativistic motion of the 
point charge, using relativistic transformations. 
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